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Abstract. We present spectra of Comet Hale-Bopp (C/1995 01) covering the range 1.4- 
2.5 fim that were recorded when the comet was 7 All from the Sun. These show broad 
absorption features at 1.5 and 2.05 fim. We show that some, but not all, of this absorption 
could be matched by an intimate mixture of water ice and a low albedo material such as 
carbon on the nucleus. However we recognise that it is more likely that the ice features 
are produced by scattering from icy grains in the coma. The absence of absorption at 1.65 
fim suggests that this ice is probably in the amorphous state. An unidentified additional 
component may be required to account for the downward slope at the longwavelength end 
of the spectrum. 



1 Introduction 


Although water is a major constituent of comets and has been detected in the gaseous phase 
a number of times (eg. see the review by Mumma et al . , 1993) the detection of solid phase 
water in comets has proven elusive. Direct observations of comet nuclei are usually prevented 
by their optically thick comae, and comets at large heliocentric distances, where the comae 
maybe thin or absent, are usually too faint for spectroscopic studies in the near-infrared 
(JHKL: 1.2-3. 6 /xm) regions where water ice is most readily detected. 

Water ice is not however expected to be confined to the surface of cometary nuclei. As gas 
sublimes from the nucleus of an active comet grains are ejected at velocities of a few hundred 
m sec -1 to 1-2 km sec -1 . If volatile icy grains are ejected in sufficient quantity they may form 
an icy grain halo that is detectable spectroscopically before the individual grains evaporate. 
The rate of evaporation of these icy grains depends upon the heliocentric distance at the 
time of ejection. Energy balance calculations (Hanner 1981) show that unless icy grains are 
very pure, they cannot survive within about 2.5 AU of the Sun, precisely the region where 
comets usually brighten sufficiently to be observed easily. Only in the case of comets which 
are bright at large heliocentric distances is ice likely to be detectable, but such opportunities 
are rare. One possible detection of ice was reported in comet Cernis (Hanner 1984) where an 
absorption feature at 2. 9-3.0 gm was deduced using CVF spectroscopy at three wavelengths, 
one in the ice band and two outside it. Campins et al. (1983) reported the detection of an 
absorption feature at 3.25 gm in comet Bowell which they attributed to a halo of icy grains 
around this comet, a possibility supported by A’Hearn et al. (1984) who reported that the 
OH production rate in Comet Bowell reached a maximum value beyond 4 AU at which time 
it was 10 times the rate recorded at perihelion. This surge in OH production was attributed 
to the sublimation of icy grains in the coma while the comet was still approaching the Sun. 



However, despite these two pieces of evidence for icy grains in Comet Bowed Campins et al 
(1982) did not detect absorptions in the 1.25-2.3/xm region which might have been expected 
from icy grains. To investigate this they investigated the scattering properties of various 
grain sizes via Mie scattering calculations and concluded that the best fit to their data was 
for particles with radii in the range 10^m and albedoes of about 5Such grains would be very 
dark for pure water ice but Clark (1981 ) has shown that even small amounts of dark material, 
eg soot, can dramatically darken icy surfaces until they have albedoes close to that of the 
dark contaminant. Campins et al (1982) thus proposed that the icy grain halo in comets 
at large heliocentic distances comprised a volatile matrix of water ice paticles in which were 
embedded dust grains which lowered the albedo to the order of 5 

Comet Hale-Bopp was discovered on 1995 July 23 ( IAUC 618 7), at which time it was at 
a heliocentric distance of 7.4 AU and a geocentric distance of 6.4 AU. It was soon realised 
that the comet was unusually bright for an object at such a large heliocentric distance. Jets 
were observed in late 1995 by a number of observers. Sekanina (1996) has discussed the 
outburst history of comet Hale-Bopp and has compared its performance to other comets at 
large heliocentric distances. Jewitt et al. (1996) and Biver et al. (1996) detected emission 
from the CO which was presumably driving the production of the coma. In this paper we 
present near infrared spectroscopy during early September 1995, when the comet was 7 AU 
from the Sun, and report the detection of shallow absorption bands at 1.5 and 2.05 /on, 
which we attribute to H 2 0 ice. 


2 Observations 

Comet Hale-Bopp was observed on five nights between UT 1995 September 4 and Septem- 
ber 9 using the cooled grating array spectrometer CGS4 of the United Kingdom Infrared 



Telescope (UKIRT). CGS4 was equipped with its 75 line per mm grating and its short focal 
length camera, giving a pixel size of approximately 1.23 by 1.23 arc sec. The slit length 
was 90 arc sec. On some nights observations were made with a 1 pixel wide slit and on 
others a 2 pixel wide slit was used, see Table 1. The spectra generally were sampled three 
times per resolution element, and the detector was stepped over two pixels to ensure full 
spectral coverage in the presence of any bad pixels. The slit was oriented N-S during the 
observations which, since the tail position angle at the time was 86 degrees, placed the slit 
almost perpendicular to the Sun vector. Object-sky pairs of observations were made with 
the telescope being nodded along the slit between observations. The amplitude of this nod 
varied between 25 rows (approx 30 arc sec) and 10 rows (approx 12 arc sec) on different 
nights. The object and sky frames were differenced and the differences co-added to form 
a final frame containing two sky-subtracted comet spectra. On-chip single exposure times 
varied from 12 to 20 seconds in the J,H and K bands and 0.4s in the L band. 

Comparison spectra were taken at airmasses that matched those of the comet and were used 
for removal of telluric features from, and flux calibration of, the comet spectrum. In the case 
of the data from September 4, for which the comparison star (CMC314447, V=9.04) was of 
spectral type A0, the final comet spectrum was divided by the mean stellar spectrum and 
the result was multiplied by a blackbody of temperature 9850K, the effective temperature of 
a star of that spectral type. We note that A0 stars have strong Brackett gamma absorption 
which, when the star is used for atmospheric cancellation, manifests itself as a spurious 
emission feature at 2.166 pm: this can be seen in the lower spectrum of Fig 1. The same 
procedure was used for the remaining nights using the calibration star BS6836 for which 
the bright star catalogue gives V=6.4 and spectral type GOV. Assuming the colours given 
by Johnson (1966) we used J =5.37 and K=5.05 for BS6836. Wavelength calibration was by 
reference to an argon emission lamp spectrum taken at the appropriate grating position each 



night. 


Details of the observations are given in Table I and the extracted spectra for each night are 
shown as Fig 1. In the case of the H band spectra, for which the spectra extracted from the 
two different nod positions showed some differences in shape, both spectra are shown (la,b) 
For the remainder we show the average of the spectra extracted from the two nod positions 
as lc-lf. 

The optics of CGS4 produce a curvature of the spectrum across the detector array which 
must be removed before extracting spectra from specific rows of the array. The curvature is 
first determined by examining the spectrum of a point source, in this case the standard star. 
The images containing both the standard star and the comet are then rebinned to remove 
the effects of spectral curvature using routines in the Starlink Figaro reduction package. In 
all the images it was clear that strong continuum emission from the comet extended over 
several rows and the nodding along the slit may not have been sufficient to clear the coma 
completely. Examination of the spatial information provided along the slit showed that 
although the flux from from the comet was assymetric, falling away more rapidly to the 
south, it was fairly sharply peaked and so in each spectrum, irrespective of the distance 
nodded, only the central five rows were extracted. This corresponds to approximately the 
full width at half maximum of the cometary signal. With the 1.23 arc second pixel size and 
a one pixel wide slit, the extracted rows are equivalent to a retangular 1.23 by 6.2 arc second 
aperture which corresponds to a projected area of 5700 by 28500km at the distance of the 
comet. In order to provide an absolute value of the flux from the comet on each date, for 
each spectrum we present the flux at the centre of the standard JHKL filter passbands which 
was obtained by extracting from the flux calibrated spectra the average flux from a spectral 
region 0.04 fim wide centered on the nominal filter central wavelength (1.25. 1.65, 2.2 and 



QQ firn respectively). These fluxes are are presented in Table 1 where it should be noted 
that not all are for the same effective aperture, since the slit width used was not the same 
for each observation, and fewer rows were extracted from the L spectrum to minimise the 
background noise. 

Throughout the observations the comet was passing through the galactic plane close to the 
galactic centre and as a consequence other stars occasionally passed through the slit. Where 
these seemed likely to contaminate the spectra, individual object-sky pairs were removed 
from the final dataset before extracting the spectra. 


3 The Composite Spectrum of Comet Hale Bopp. 

Since the spectral features of interest are broad, and lie close to the edge of atmospheric 
windows, various wavelength settings were used in order to define the adjacent regions of 
continuum as accurately as possible. To improve the overall signal to noise ratio, all the 
spectra were combined to produce a single composite spectrum. The original flux calibrated 
spectra were divided by the incident solar flux which was calculated from Smith and Gottlieb 
(1974). The three K band spectra taken on September 4, 5 and 9 were scaled to unity at a 
wavelength of 2.2 /xm. The H band spectra were then scaled to the mean K band spectrum 
at 1.66 [im. Data points with negative or zero flux values, believed to be due to noisy pixels 
and cosmic ray events, were deleted from the files as physically unrealistic. 

All the data were then resampled to a common set of wavelength values spanning the entire 
range observed. This was done using a Gaussian convolution with a FWHM of 0.0009 p,m. 
This corresponds to 1.5 spectral channels for the nights with the highest spectral resolution 
and results in a slight smoothing of the overall spectrum. A weighted average was calculated 



from the convolved results. In order to determine the weighting we evaluated the point- 
to-point scatter in the scaled reflectance values for each observation over a 0.1 /xm range 
which best represented a neutral part of the spectrum. Regions of high telluric absorption 
and the region of the 1.5 and 2.05 /xm absorptions were excluded. Based upon the standard 
deviation of the point-to-point scatter the spectra were given the weights listed in Table I. 
After calculating the average spectrum, data points with standard deviations greater than 
0.15 were deleted from the spectrum. 

The J band spectrum was featureless and had no wavelength overlap with the remainder of 
the spectra taken. It is not plotted here. The L band spectrum (also not plotted) was very 
weak and even after combining several rows and smoothing the result revealed only a vague 
impression of a broad absorption between 3.0 and 3.2 /xm with some hint of continuum from 
3. 2-3. 6 /xm. From this we estimate the L band continuum flux density to have been of the 
order 1 x 10 -16 Wm -2 in an aperture 1.2 by 4.8 arc sec oriented EW x NS. 

The individual H and K spectra are shown in Fig 1 and the resulting composite 1.4-2. 5 /xm 
spectrum is shown in Fig 2. In both figures two broad absorption features can be seen centred 
at wavelengths of 1.5 and 2.05 /xm. In addition the spectral reflectance drops steeply at the 
long wavelength edge. Believing these absorptions to be due to the presence of water ice, we 
have investigated by means of Hapke theory, whether they can be explained by material on 
the nucleus of the comet. Because of limitations in our modelling code, we do not address 
the issue of an icy grain halo in this paper although we recognise that the scattering by 
grains must contribute significantly to the observed flux since the reflected light at 2.2 /xm 
from a 10km nucleus of albedo 0.04 would be only 8 x 10 -19 Wm -2 /xm -1 , only 0.1% of the 


flux observed. 



4 Mathematical Modeling - Equations and Assump- 


tions 

4.1 Intimate Mixtures 

Hapke (1981, 1986, 1993a, 1993b) has developed a series of equations, based on a two-stream 
approximation to the equations of radiative transfer. This approach provides the ability to 
calculate the reflectance of a particulate surface given assumptions or knowledge about its 
physical nature and the solid phases comprising it. Such mixtures are commonly referred to 
as intimate mixtures where individual components exist as separate grains that are in turn 
thoroughly mixed at the micrometer scale. 

The bidirectional reflectance of a surface, at each wavelength (A), can be expressed as 

r(\J,e,g,w,h,S(0)) = ^-^-{[l + B(g)}P(g) + H(ti)H(ti o )-l} (1) 

where i and e are the angles of incident and emergent light, respectively (fi 0 = cos i and 
/i = cos e), g is the phase angle between i and e, w is the average single scattering albedo of 
the particles composing the surface, h and S(0) are related to the width and amplitude of 
B(g) (the opposition surge) that is due to the microstructure of the surface (e.g., porosity, 
particle size distribution, and compaction rate with depth), P(g) is the the average particle 
phase function, and H(/x 0 ) and H(/z) are Chandrasekhar’s (1960) H-functions. In all the 
models discussed below we assume an opposition effect, that is is similar to that of other icy 
surfaces (i.e. h=0.05), an isotropic phase function (i.e. P(g)=l), and no internal scattering 
within the particles. The average single scattering albedo is related to the particles’ relative 
mass fractions, solid densities, and wavelength dependent volume absorption coefficients (a) 



(Hapke 1981, equation 17). The absorption coefficient is related to the imaginary index of 
refraction ( k ) through the dispersion relation (a = Airk/ A). See Roush et al. (1990, 1996), 
Roush (1994), and Cruikshank et al. (1996) for further details regarding the modeling. 

Because absolute albedo information cannot be deduced from the observational data, due 
to the lack of knowledge of the diameter of the comet, both the observational data and the 
calculated spectra of the resulting mixtures are normalized by scaling them to a value of 1.0 
at a wavelength of 1.73 fim for comparison purposes. 

4.2 Spatial Mixtures 

Spatial mixtures represent a surface having distinctive albedo regions present on a scale 
smaller than the spatial resolution of the observations. These are modeled by a linear 
combination of the reflectances of each region weighted by its spatial extent. The total 
reflectance of the surface ( Rj ) is expressed mathematically as 

Rj = A* Ra + 1 B*Rg + ■ • • .V * Rpj ( 2 ) 

where A. B, and N are the spatial coverage and Ra , Rb , and Rn represent the spatial 
coverage and reflectance (as calculated from equation 1) of each component, respectively. 

4.3 Optical Constants 

The optical constants of water ice and amorphous carbon were used to calculate the re- 
flectances discussed below. Water ice was selected because it exhibits absorption features at 
the same wavelengths as the observational data and because it is expected to be present in 
comets. Amorphous carbon was selected because it is a low albedo, spectrally neutral ma- 
terial. The real and imaginary indices of refraction of hexagonal water ice at LOOK (Roush. 



1996) and amorphous carbon (Rouleau and Martin. 1991) were used since this corresponds 
closely to the expected black body temperature at the distance of the comet. These refrac- 
tive indices were linearly convolved to the same spectral resolution as the observational data. 
As mentioned above, the reflectances are directly influenced by any changes associated with 
the optical constants used in a calculation. If the optical constants used in these calculations 
are revised or supplanted by newer more accurate values then the results and conclusions 
presented below will be affected accordingly. 


5 Calculations and Comparisons to the Observational 
Data 

5.1 Pure Ice 

Fig 3 shows the resultant calculations for a variety of grain sizes of pure water ice and a 
comparison of these to the comet Hale-Bopp data. The grain sizes cover a range that produce 
absorption features spanning the strengths observed for comet Hale-Bopp. However, it should 
be noted that the two finer grain sizes may violate the assumption of Hapke theorv that the 
particles’ cross sectional areas are greater than the wavelength of the incident radiation 
(xd/\ >> 1). The H 2 0 ice spectra in Fig 3 exhibit absorption features centered near 1.5- 
and 2.0 /cm and the wing of a strong feature centered near 3 fim. It is readily apparent 
that a surface composed solely of water ice can roughly reproduce the water ice features 
seen in the comet Hale-Bopp spectrum at the shorter wavelengths but cannot reproduce the 
relative strength of the three bands. In particular the 3 fim band is too strong relative to 
that expected from the model spectra. This suggests that some additional component or 
alternative mechanism gives rise to the extra absorption at wavelengths >2.2 [im. 



5.2 Intimate Mixtures 


Because water ice cannot simultaneously reproduce the relative strengths of the 1.5, 2.0, 
and wing of the 3 /im absorptions, only data lying between 1.4 and 2.2 //m were used to 
model the observational data. The model represented the surface as an intimate mixture of 
water ice and amorphous carbon and the fitting algorithm consisted of a downhill simplex 
approach (Press et al. 1989) that minimized the expression. 


x 2 = E 



(3) 


where the sum is over wavelength, r 0 and r c are the observed and calculated reflectances, 
respectively, and a 0 is the standard deviation of the observed reflectance. This approach was 
used to determine the relative mass fraction and grain size for each component. Once the 
algorithm converged to a particular solution, the resulting parameters were used to calculate 
the scaled reflectance over the entire wavelength range of the observational data. In an 
attempt to avoid local mimima in the parameter space, the algorithm was initialized at 
several different locations; only two of the best solutions are illustrated in Fig 4. 


Fig 4a shows the results of the spectrum fitting procedure for a homogenous surface com- 
posed of roughly one quarter water ice and three quarters amorphous carbon having grain 
diameters of 133 and 14 /im, respectively. The relative strengths of the 1.5 and 2.0 \im 
bands are reproduced well by the model, however, the absorption at longer wavelengths is 
not. We note a discrepancy between the model and the comet spectrum at 1.65 //m. At this 
wavelength a temperature-sensitive absorption band of crystalline ice occurs. The absorp- 
tion is minimal for T > 200K, and increases with decreasing temperature (Fink and Larson 
1975). Amorphous ice at T=100K shows minimal absorption at this wavelength (Fink and 
Sill 1982). The model spectrum presented here shows absorption at 1.65 f.im because the 



optical constants used were derived from spectra of crystalline ice at 100K in which moder- 
ately strong absorption occurs at that wavelength. The absence of the absorption feature in 
the comet spectrum is therefore suggestive of the presence of amorphous ice. 

The exothermic phase transition from amorphous to crystalline water ice occurs in the tem- 
perature range 125-150K (Mumma et al. 1993) and so amorphous ice might not be expected 
close to the surface of comets which have made one or more perihelion passages at distances 
less then a few AU. However, while comet Hale-Bopp is not new in the Oort sense, the 
detection of CO shows that some very volatile materials have been retained in the nucleus 
and are only now becoming manifest. Amorphous or cubic phases of water ice could also 
have survived and may now be being brought from within the comet by the jets observed 
in the visible. At present optical constants of these phases of water ice are not available 
at the wavelengths associated with the comet data presented here but when such data are 
available it may be that models using such values might provide a better correspondence to 
the observations. The alternative possibility that the ice is crystalline but at T > 200K is 
highly improbable for a comet at 7 AU from the sun. 

A particle size distribution was investigated by including an additional water ice component 
having a different grain size than the initial water ice component. A resultant solution is 
shown in Fig 4b. In this mixture, water ice occurs as two grain sizes, 133 \im and 27.7 mm 
that combined have greater abundance («40%) while the amorphous carbon has a larger 
grain size 21 (im. but is less abundant than in the two component mixture. However, the 
spectral characteristics of the water ice with grain sizes of =s28 mm would closely resemble 
those of amorphous carbon, i.e. it would have a low albedo and be spectrally neutral over the 
entire wavelength region investigated here. Nevertheless, the resulting mixture reproduces 
the 1.5 and 2.0 fim band strengths well, but again cannot produce the long wavelength 



absorption. 


5.3 Spatial Mixtures 

The ability of spatial mixtures of water ice and amorphous carbon to describe the observa- 
tions were also investigated. Because amorphous carbon is strongly absorbing its calculated 
reflectance is independent of grain size (at least over the grain size range that remains consis- 
tent with the assumptions of Hapke scattering models). The reflectances were calculated for 
pure water ice having grain diameters of 5-, 10-, 50-, 100-, 500-, and 1000 fmi. A commercially 
available curve fitting routine ( SigmaPlot ) was used to determine the relative abundances of 
mixtures of amorphous carbon and water ice. This software relies on a Marquardt-Levenberg 
algorithm. Like the intimate mixtures discussed above, only the observational data lying be- 
tween 1.44 and about 2.2 fim that were scaled to 1 at 1.73 fim were used in the fitting 
procedure. Once a solution was determined, the derived relative abundances were used to 
calculate a spectrum over the entire wavelength range of the observational data and equation 
3 was used to calculate \ 2 . 

The best fitting spatial mixtures, shown in Figs 4c and 4d, resulted from the two finest 
grain diameters of water ice and yielded similar relative abundances. Both consist of 99% 
amorphous carbon and 1% water ice. The \ 2 values were very similar to the best fit intimate 
mixtures. 


6 Conclusions 


Due to limitations in out modelling code we have restricted this analysis to intimate and 
spatial mixtures on a solid surface although we recognise that in reality reflection from the 



nucleus is not the major contributor to the flux recieved from the comet even at this large 
heliocentic distance. Least-squares fitting proceedures that model the reflection of incident 
sunlight by a particulate surface suggest that water ice can reproduce the relative strengths 
of the 1.5 and 2.0 /.im absorptions observed in comet Hale-Bopp's spectrum. A range of 
solutions, involving one or two particle sizes for water ice. intimately mixed with a low 
albedo neutral color material such as amorphous carbon, suggest that the grain size of water 
ice responsible for the observed band strengths is ^100 firri and constitutes a relative mass 
fraction of roughly 15-25%. Spatial mixtures of amorphous carbon and water ice favour a 
surface dominated by amorphous carbon (99%) with minor areas of fine grained water ice 
(1%). Spatial mixtures and intimate mixtures provide solutions that, based on the quality of 
fit to the observational data, are comparable. However, none of the mixtures can reproduce 
the observed absorption at wavelengths >2.2 /im. This suggests that either (1) another 
material is present that gives rise to the additional absorption but lacks other distinctive 
spectral features at these wavelengths; or (2) some other mechanism is giving rise to the 
additional absorption, e.g., absorption due to icy grains present in the coma of comet Hale 
Bopp. Scattering from particles in an icy grain halo may account for the discrepancy at 
the longer wavelengths. The presence of a large number of neutral scattering grains in the 
coma may also have the effect of diluting the absorption features seen in a relatively large 
beam. This would reduce the amount of neutral material required in each grain to match 
the relatively shallow features seen in the data. We are presently working on modeling code 
to explore scattering models. Although optical constants at the relevant wavelengths for 
amorphous ice are not presently available, we are attempting to obtain such data in order to 
explore the possibility that the absorptions in comet Hale Bopp are due to ice in the vitreous 
phase. 
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fim 

R 

\/S\ 

Sec z 

Nod 

arcsec 

Aperture 

arcsec 

Flux 
at A 

fim 

Flux 

Wm ~ 2 fim~ x 

sp 

wt 

4.3 

2 

1.85-2.51 

450 

1.65 

30.7 

2.48 * 6.2 

2.2 

3.39 ± 0.10 x 10“ 15 

1 

4.3 

2 

1.32-1.98 

320 

1.76 

30.7 

2.48 * 6.2 

1.65 

3.78 ± 0.29 x 10~ 15 

1 

5.2 

1 

1.84-2.50 

850 

1.57 

12.3 

1.24 * 6.2 

2.2 

2.19 ± 0.16 x 10~ 15 

1 3 

6.2 

1 

1.01-1.36 

450 

1.58 

12.3 

1.24 * 6.2 

1.23 

8.82 ± 0.31 x 10 -15 

0 

7.2 

1 

1.57-2.23 

750 

1.58 

12.3 

1.24 * 6.2 

2.2 

2.23 ±0.19 x 10~ 15 

3 

9.2 

1 

1.64-2.30 

770 

1.64 

18.4 

1.24 * 6.2 

2.2 

2.64 ± 0.29 x 10" 15 

5 

9.3 

1 

3.01-3.69 

1300 

1.56 

18.4 

1.24 * 4.9 

QQ 

3.0 ± 1.0 x 10~ 16 

0 


Table 1 : Details of the Observations taken 

8 Figure Captions 

Figure 1 The six separate observations used to determine comet Hale-Bopp’s composite 
spectrum. The spectra in a and b were scaled to the mean K-band spectrum at 1.66 fim. 
The spectra in c-f were scaled to 1.0 at 2.2 /xra. 

Figure 2. The composite spectrum of Hale-Bopp (top curve) and associated standard devi- 
ations (bottom curve). The data have been resampled to common wavelength scale, points 
with standard deviations greater than 0.15 have been deleted, and rescaled to 1.0 at 1.73 
//m. 

Figure 3. Comet Hale-Bopp data (heavy solid line) compared to the reflectance calculated 
for water ice having grain diameters of 1 fim (dashed line), 5 fim (dotted line), and lOfirn 
(thin solid line). 




Figure 4. Comet Hale-Bopp data (solid lines) compared to (a) the reflectance calculated 
for a two component intimate mixture of water ice and amorphous carbon, (b) a three 
component intimate mixture containing two particle sizes of water ice. (c) a two component 
spatial mixture of amorphous carbon and 5 /jm diameter water ice, and (d) a two component 
spatial mixture of amorphous carbon and 10 yum diameter water ice. 
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Reflectance Scaled to 1 .0 at 1 .73 ^im 
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